Kondo time scales for quantum dots — response to pulsed bias potentials 
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The response of a quantum dot in the Kondo regime to 
rectangular pulsed bias potentials of various strengths and du- 
rations is studied theoretically. It is found that the rise time is 
faster than the fall time, and also faster than time scales nor- 
mally associated with the Kondo problem. For larger values 
of the pulsed bias, one can induce dramatic oscillations in the 
induced current with a frequency approximating the splitting 
between the Kondo peaks that would be present in steady 
state. The effect persists in the total charge transported per 
pulse, which should facilitate the experimental observation of 
the phenomenon. 

PACS numbers: 72.15.Qm, 85.30.Vw, 73.50.Mx 



The theoretical predictions p]] of consequences of the 
Kondo effect for the steady state conduction through 
quantum dots began a decade ago. At low temperatures, 
a narrow resonance in the dot density of states can form 
at the Fermi level, leading to a large enhancement of 
the dot's conductance, which is strongly dependent on 
temperature, bias, and magnetic field. Many of these ef- 
fects have been recently observed by a set of beautiful 
experiments by several groups [||. These successes, sup- 
plemented by the anticipation that time dependent ex- 
periments H are not far behind, have spurred a number 
of theoretical groups M to consider the effects expected 
when sinusoidal biases or gate potentials are applied. 
Surprisingly, the application of steps or pulses, which can 
provide a less ambiguous measure of time scales, have not 
been considered until very recently [jsj; the latter work 
considered the time dependent change in linear response 
conductance when a stepped potential was applied to a 
gate, thereby shifting the dot into the Kondo regime. In 
the present work we consider the response of a dot al- 
ready in the Kondo regime to a sudden change of the 
bias potential across the dot. We show that the physics 
is qualitatively different from the latter case leading 
to a different range of characteristic times and physical 
phenomena. 

While Ref. H studied the time scale for the system 
to go from one equilibrium configuration to another, we 
study here the time scale to go from an equilibrium con- 
figuration to a non-equilibrium one; and then back to 
equilibrium again. We find that these latter two times 
scales are very different from each other, the first being 
much shorter - and also much shorter than the time scale 
of [B. Furthermore, if one applies a rectangular bias 



pulse large enough to split the Kondo resonance, then 
there appear current oscillations || at a frequency char- 
acterizing the splitting between the Kondo peaks. We 
show that these current oscillations cause oscillations in 
the charge transported through the dot as a function of 
pulse length, thus providing a clear experimental signa- 
ture. The damping of the oscillations provides an addi- 
tional time scale which also can be measured directly. 
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FIG. 1. The instantaneous current I(t, r) induced by a bias 
pulse of height evo =0.01 beginning at t — and ending at 
t = t — 600 (heavy solid line). The portion for t > r is 
translated back to the origin and inverted (light solid line) 
so that the rise time and fall time can be easily compared. 
The area between these curves represents the excess charge 
forced through the dot above what would flow [G(vo)t] if 
the response to no were instantaneous. The value of vq is 
small, and these curves represent the very beginning of the 
non-linear response regime. As vo increases further, the rise 
time becomes shorter [oc 1/vo], while the fall time remains 
long. The dashed line represents the extent to which the 
Kondo state would be formed (at zero bias) if the Kondo 
coupling were to be suddenly turned on at t = (see text). 
Inset: Schematic of the quantum dot during a pulse. 

We model the quantum dot by a single spin degenerate 
level of energy edot coupled to leads through tunnel bar- 
riers, as illustrated schematically in the inset to Fig. [j. 
The Coulomb charging energy U prevents the level from 
being doubly occupied. We apply a time-dependent po- 
tential across the leads v(t) = vq for < t < t and 
zero otherwise, so that the energy Ek of an electron in a 



lead is time dependent. Here k represents all lead quan- 
tum numbers other than spin, including the labels, left 
and right, specifying the lead in question. Then we have 
Ek = E k (t) = e k ± hev{t), the sign depending on which 
lead k refers to, where e k is the band energy. The system 
may be described by the following Anderson hamiltonian: 
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with the constraint that the occupation of the dot cannot 
exceed one electron. Here cj. creates an electron of spin a 
in the quantum dot, with n CT the corresponding number 
operator; c' k<T creates an electron in the leads. 

The general features of the static equilibrium spec- 
tral density when the dot level e is sufficiently below 
the Fermi level (taken at zero here) are well known. 
There is a broad resonance of half-width ~ Tdot at 
an energy ~ edot plus a sharp temperature sensitive reso- 
nance at the Fermi level (the Kondo peak) , characterized 
by the low energy scale Tr (the Kondo temperature), 
T K ~ Z?'exp(-7r|e do t|/r dot ), where D' is a high energy 
cutoff ||. Throughout this work energies and temper- 
atures are given in units of Tdot, and times in units of 
1/Tdot, with h = 1. Specifically, the calculations in this 
Letter were made with edot = —2, which leads to a Kondo 
temperature Tk ~ 0.0025. We present results only for a 
dot with left-right symmetry. 

The current into the dot depends on the time t and 
paramctrically on the pulse length r, and is given by 
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It may be divided into contributions Ii e tt(t, r) and 
bright (^jt) by respectively restricting the k summations 
to the appropriate lead. The transport current is 
then I{t,r) = |[Ji e ft(i,r) — /right (t, r)]. We calculate 
the Keldysh propagators corresponding to the angular- 
bracketed expectation values in (||) for each lead, and 
hence obtain I(t,r). A more experimentally acces- 
sible quantity, the total charge transported, Q(t) — 
J °° dt I{t, t) is also calculated. Our calculations use the 
non-crossing approximation (NCA), which is reliable for 
temperatures down to T < Tk ||. Its time-dependent 
formulation has been developed elsewhere |[o| . The finite 
bias on the leads is taken into account by introducing a 
time-dependent phase in V k in (^). 

In Fig. |l|, we show the response to a long, weak bias 
pulse, which illustrates some features that will persist 
for still stronger pulses. The most important feature is 
that the rise time of the current is shorter than the decay 
time, leading to a greater charge transport than would 
have occurred if the response at each end of the pulse 
had been instantaneous. 

Another comparison can also be made: how does this 
time scale for turning on the current compare with how 



long it takes to form the Kondo state itself? The latter 
was obtained in Ref. || , by calculating the turn-on of the 
linear response current when the parameter edot was sud- 
denly switched into the Kondo regime. This is illustrated 
by the dashed curve in Fig. |l|, which is the calculation 
of Ref. H for the present parameters, with the abscissa 
scaled to make the curves coincide at saturation (large r 
and large t, with t < r). It is clear that the formation 
time for the Kondo state is significantly longer. 

A closer look at Fig. |l] reveals that approximately a 
third of the initial rise (in both curves) when the bias is 
turned on (and the same fraction of the initial fall when 
the bias is turned off) occurs essentially instantaneously 
on the scale of the graph. This fast response occurs on the 
time scale to = 1/Tdot- This represents the part of the 
conductance arising from the fact that the dot level has 
finite width which overlaps the Fermi level, and would be 
present even in the absence of the Kondo peak. 

The Kondo derived features from switched pulses can 
be partially understood by reference to the steady state 
spectral densities for a dot with a fixed bias equal to vq. 
Once eva ^> max(T^, T), the Kondo peak becomes split, 
with a peak separation associated with the potential 
±t}Vo of each respective lead, as discussed by Wingreen 
and Meir . Such spectral densities for increasing values 
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FIG. 2. Steady state spectral densities for different values 
of the bias. The left panel shows the whole spectral function 
for the bias equal to the vo of Fig. [jj with the Kondo region 
expanded in the inset; the equilibrium spectral function is also 
shown. The three right panels show just the Kondo region 
for increasing biases. The heavy curves represent the total 
spectral density, while the light curves represent the occupied 
portion. The oscillations that we find in the next two figures 
presumably represent transitions from the lower Kondo peak 
in the occupied spectra to the upper peak in the unoccupied 
spectra (i.e., the difference between the pairs of curves). Note 
that the vertical scale factors for all four panels (but, not the 
inset) are identical. T = 0.0015 everywhere. 
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vq are illustrated in Fig. 0. This splitting should provide 
a sharp excitation frequency u ~ evo for the dot, with 
an electron being excited from the lower peak at energy 
— \v to the upper peak at + \v$. We should therefore 
expect current oscillations of this frequency to be induced 
by the rapid turn on of the voltage pulse. For larger vq's 
than in Fig. ^ the frequency u provides a fast charac- 
teristic rise time t\ = tt/2ui, and we may expect that the 
rise time depicted in Fig. 0, is a precursor to this effect. 
Fig. |3|, showing the results for a much larger vq, not only 
confirms this interpretation, but shows how dramatic the 
current oscillations can be. In what follows, we refer to 
them as split Kondo peak (SKP) oscillations. 

How can one experimentally distinguish SKP oscilla- 
tions from those reported jnj for the resonant level model 
(U = 0), which occur as each lead adjusts to its new 
chemical potential (NCP)? Our calculations show that, 
unlike the SKP oscillations which occur robustly at a 
frequency ~e«o, the NCP oscillations are sensitive to the 
details of the exciting pulse. First, with the resonant level 
at the Fermi level to simulate a Kondo peak, the NCP 
oscillation frequency is equal to evo only in the special 
case that a pulse of magnitude |et>o| is applied to a single 
lead, with the other held at the original Fermi level. For 
equal amplitude pulses of ±evo/2 to the respective leads, 
the oscillation frequency is, on the other hand, evo/2. 
In a more general case of asymmetrically applied pulses 
there are two oscillation frequencies, one associated with 
the new chemical potential of each lead. If in addition, 
the resonant level is not at the original Fermi level, the 
oscillations are even more complicated, with multiple fre- 
quency components. 

The NCP oscillations are damped on the time scale 
~ 1 /Tdot- When this becomes shorter than the oscillation 
period, the oscillations become overcritically damped 
and do not appear. We find that they only occur for 
e^o/Fdot > 1, and not for the much smaller values of this 
ratio used in the large U examples shown in this paper. 
SKP oscillations, on the other hand, occur only when 
e^o/Tdot *C 1, and hence not simultaneously with NCP 
oscillations. This explains the robustness of the oscilla- 
tions in the Kondo regime. 

SKP oscillations can probably be excited by other 
types of applied waveforms, and a parametrically driven 
form of them is probably responsible for the resonances 
reported by Schiller and Hershfield H. However, we sus- 
pect that the well-defined robust frequency for them pro- 
vided by a square pulse, plus the lack of one for NCP 
oscillations, may provide the key for unambiguous ex- 
perimental identification. 

Fig. ^ shows an additional feature — the damping of 
the SKP current oscillations. This damping presumably 
arises from the fact that the split Kondo peaks in the 
spectral density have widths, giving rise to dephasing. 
These widths are substantially greater than the width of 
the Kondo peak in equilibrium, as the application of a 
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FIG. 3. The instantaneous current I(t,r) induced by a 
strong pulse (an order of magnitude larger than the width 
of the zero bias Kondo peak) for several pulse lengths r. The 
units of both r and t are both 2ir/evo, which is almost pre- 
cisely the "period" of the induced oscillations (although we 
see the frequency slightly decrease from uj — evo as the bias 
is lowered towards Tk). Note the reversal, for certain pulse 
lengths, of the direction of the instantaneous current after the 
end of the pulse. 

bias begins to weaken the Kondo effect. The reason for 
this has been pin-pointed by Wingreen and Meir p| , in a 
model perturbation theoretic calculation whose essential 
feature is that the particle-hole phase space restriction 
from the Pauli principle in the leads for incoherent spin 
processes is eased. For biases evo » T, this gives a rate 
enhancement factor of ~ cvq/T relative to the zero bias 
case. For T < Tk, one should replace T by Tk in the 
above argument; in this case the width of the Kondo 
peak(s) would be proportional to evo for evo 3> Tk, but 
would saturate at ~ Tk for evo <C Tk- Translating this 
into the damping time t^ of the oscillations, gives ti cx 
1/uo at large biases, but with a saturation t% ~ 1/Tk 
when evo < Tk- The expected saturation as the split 
Kondo peaks merge removes the oscillations for evo < Tk 
(cf. FigJ). 

Fig. [|also shows the fall off in current after the voltage 
pulse has been turned off. After an initial drop due to 
the fast non-Kondo time scale t = 1/Tdot, the current 
falls very slowly. Its decay does not follow single time 
scale, but it is clear that the total time required for the 
decay is very long. While the initial fall-off after to is bias 
dependent and relates to the time t\, the long tail has a 
characteristic time t% ~ l/?k, the longest time in the 
problem. This is to be expected: the system is trying to 
regain an equilibrium where the Kondo peak is no longer 
split or broadened. 

Although we only show the results for one tempera- 
ture, T — 0.0015 ~ Tk, we have made calculations for a 
number of higher and lower temperatures. The results 
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FIG. 4. 1/vo times the derivative with respect to pulse 
length r of the total charge Q transported during and af- 
ter a pulse vs. r. The curves go to the zero bias conductance 
G(0) at r = 0, and each go to the respective values of the 
finite bias conductance G(vo) for large r. 

are not very sensitive to temperature when T < Tk- At 
higher temperatures, the oscillation amplitude decreases, 
and each time ti, ti, or £3 shortens when T gets large 
enough to equal the reciprocal of the time in question 
(T ~ 1/ti). The fall-off time and oscillation decay time 
t-i are the first ones affected. At T ~ evo, the rise time t± 
also shortens. However, the nonlinear effects persist to 
temperatures much higher than Tk presumably because 
of the logarithmic nature of the Kondo peak. 

How can these effects be seen experimentally? Wc 
suggest that the quantity Q(t) (the charge transported 
through the dot by a single pulse of length r) will con- 
tain sufficient memory of the charge already transported 
at the time r when the bias is turned off, to directly re- 
flect the current oscillations and their damping. Q(t) 
is proportional to the average current transported by a 
series of successive pulses, provided that the time be- 
tween pulses is 3> ta; it thus corresponds to an easily 
measured quantity. To restore the features to the promi- 
nence which was there before the time integration, it is 
suggested that dQ(r)/d,T be measured. In Fig. [I] we dis- 
play our predicted dQ(r) / dr plotted vs t for a number of 
values of bias amplitude vq. The information about the 
frequency and damping of the oscillations, and hence of 
the time scales ii and t% is preserved. Similar oscillatory 
structure for fixed r occurs in plots of dQ/dvo vs vq- 

In conclusion, we have studied the response of a quan- 
tum dot in the Kondo regime to a large pulsed bias volt- 
age across the leads. We find that the rise time of the 
instantaneous current is related to the period of current 
oscillations that are set up. These oscillations have a fre- 



quency corresponding to the energy difference between 
the split Kondo peaks. The damping rate of these oscil- 
lations is related to the widths of the split Kondo peaks. 
The fall-off of the current when the pulse is turned off 
can take much longer than the rise, with a tail of length 
approaching 1/Tk- These effects may be studied ex- 
perimentally by measuring the total charge transported 
through the dot during and after a voltage pulse. 
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